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Abstract 

We study the Nambu-Jona-Lasinio model at finite temperature and finite density by using the 
functional renormalization group. The RG flows of the four-Fermi coupling constant in the NJL 
model are investigated. We obtain the chiral phase boundary in cases of the large-Y leading 
approximation and an improved approximation. The large-Y nonleading term at the vanishing 
temperature has a singularity at the Fermi surface. We show that the quantum corrections by the 
large-Y nonleading term largely influence the phase boundary at the low temperature and high 
density region. 
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I. INTRODUCTION 


One of the important subjects in elementary particle physics is to understand the basic 
properties of hadronic matter described by Quantum Chromodynamics (QCD). The phase 
diagram and equation of state of QCD have been investigated both theoretically and exper¬ 
imentally. A number of studies indicate that QCD at hnite temperature and hnite density 
has various phases such as the quark-gluon plasma phase with effective chiral symmetry and 
the color superconductivity phase. 

Investigating QCD vacuum need non-perturbative method due to strong interactions of 
QCD. Various non-perturbative methods are employed in order to study the phase structure 
of QCD. Lattice gauge simulation is very powerful because it is the hrst principle calculation 
respecting the gauge symmetry. However it is difficult to maintain the chiral symmetry on 
the lattice. The even more difficult problem which is known as the sign problem occurs at 
hnite chemical potential because the Dirac operator becomes complex. 

Other non-perturbative methods, such as the mean-held approximation (MFA), the 
Schwinger-Dyson equations (SDE) or the large-V expansion, have been applied to QCD. 
These methods do not have the difficulty of maintaining the chiral symmetry nor the sign 
problem. However, their approximations break the gauge symmetry and its systematic im¬ 
provement is not easy. 


by using the 
6| . The main 


In this paper, we study the dynamical chiral symmetry breaking (Dy;Sl^ 
functional renormalization group (FRG), another non-perturbative method 
idea of the FRG is solving the effective action Fa generated by the integration of huctuations 
with higher momentum modes A < |p| < Aq where A is the infrared (IR) cut-off scale and 
the bare action is dehned at the initial scale Aq. The evolution of the effective action with 
decreasing the IR cut-off scale is described by the functional differential equation, so-called 

irn, 


the Wetterich equation 


aFA = :rSTr 



A t 


) 




( 1 ) 


This equation itself is exact. However, we cannot solve it exactly. Therefore we solve the 
effective action in the truncated theory space. The improvement of approximation can be 
done by enlarging the theory space, which are easier and more systematic than the case of 
the MFA or the SDE. 
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The D^SB has been inves tigate d by using low energy effective models such as the Nambu- 


Jona-Lasinio (NJL) mode 


bosonization methods 


16 


15| which is described by the four-Fermi interactions. Re- 


18[ are often applied in order to obtain the physical quantities in 


macro scale from the interactions in micro scale via the FRG. In such analyses, however, the 
theory space is widely expanded and the analysis becomes highly complicated. Recently, a 
new method called the “weak solution” is proposed in order to evaluate the macro physical 
values without re-bosonization. So far, the weak solution method can be applied to a 
specihc type of the approximated FRG equation corresponding to the Fermionic large-iV 
leading approximation. 

In this paper we analyze the NJL model at hnite temperature and hnite density without 
re-bosonization method nor weak solution method. Our theory space and quantum cor¬ 
rections included can not be treated by the weak solution method up to the present. We 
compare the behaviors of the RG flows of the four-Fermi coupling constant between the 
large-leading approximation and that beyond it. We actually evaluate the RG flows of 
the inverse four-Fermi coupling constant although its physical and mathematical ground is 
not completely clear. In this manner we dehne macro physics beyond the MFA in the NJL 
model. However our method does not allows us to evaluate a critical endpoint in the phase 
diagram. Our work aims to motivate future studies of phase diagrams including the first 
order phase transition. 

This paper is organized as follows: In section [TTl we briefly explain the NJL model in view 
of FRG and our analysis method. The RG equations of the four Fermi coupling constant 
are numerically analyzed at hnite temperature and density in section IIIIl Summary and 
Discussion are given in section IIVI 


II. NAMBU-JONA-LASINIO MODEL IN VIEW OF FRG 

In this paper, we study the NJL model as a low energy effective model of QGD. The NJL 
Lagrangian is given by 

where G is the four-Fermi coupling constant and N is the number of degrees of freedom 
for the Dirac Fermion 'ip, pj. The four-Fermi interactions are generated by the QGD gauge 
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interactions. This Lagrangian is invariant under the chiral U(l) transformation: 'ip —>■ 

In present paper, we analyze the U(1 )l x U(l)i{ symmetric system to study basic structures 
and behaivors of the four-Fermi interaction at finite temperature and density. The four-Fermi 
coupling constant in the effective action corresponds to the chiral fluctuation: G {{'ip'ipf). 
Therefore, the divergence of the four-Fermi coupling constant in the course of RG flow means 
the signal of the DySB as the second order phase transition. 

The RG flow equation of the four-Fermi coupling constant for the Ad sharp cut-off scheme 
in the large-iV leading approximation is given by 

dtg = -2g + 2/ (3) 


where t is the dimensionless RG scale defined hy t = log(Ao/A), and g = GA^/dvr^ is the 
dimensionless four-Fermi coupling constant 19|. This RG equation has an ultraviolet fixed 
point g* = 1. When we solve this equation with an initial value go > 1, the RG flow diverges 
at the following finite scale. 


A = ^ log 


9o 


(4) 


.90 - 1 

We cannot solve the RG flow after tc because of this divergence. Recently, the “weak 


solution” method 


20[ | is introduced in order to define the flows after this divergence of g. 


In this section, we propose another analysis method to evaluate the RG flow of g. 

We define the RG flow equation of the inverse four-Fermi coupling constant, g = 1/g, 
and we have the RG flow equation of g as follows. 


dtg = 2g-2. 


(5) 


The solution of this equation is 


9(f) = (9„ - l)e“ + 1. 


( 6 ) 


The flow with an initial value < 1 reaches zero at R, g{tc) = 0, which corresponds to 
the divergence of g. In order to interpret the RG flow with negative g values, the Euclidean 
NJL action is bosonized. 


5'njl = / 


p){^ -h 7oA)^ ~ ^ -hp){a \ i757r)V’ 


(7) 
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where h is the Yukawa coupling constant and, a and vr are auxiliary fields. The inverse 
four-Fermi coupling constant corresponds to the mass of the mesonic potential: 


m 


2 


G 


( 8 ) 


Therefore the negative RG flow of g might indicate that the curvature of the mesonic poten¬ 
tial at the origin becomes negative and the potential takes a double-well shape, thus, DySB 
develops. Such a crude treatment has turned out to be right in case when the weak solution 
can be defined j^|. This simple analysis method follows the RG flow after DySB and reach 
the IR limit A —)■ 0. In the next section, this method is applied to the NJL model at finite 
temperature and finite density. 


III. THE RG FLOW AT FINITE TEMPERATURE AND DENSITY 


A. RG flow equations of the four-Fermi coupling constant 


We analyze the Euclidean effective action in the local potential approximation 2l|, 

G 


= J dr J d^x 




(9) 


where fd is the inverse temperature 1/T and p is the chemical potential. The evolution of 
Ta is described by the Wetterich equation. In this case it reads. 


= -Tr 


dtRA 


di,i) 


+ Ra 


( 10 ) 


where dt = — A^ and Ra is the cut-off profile function for the momentum of Fermion 
defined in Appendix]^ We introduce the simple notation, 

t , 


p(i)i) __ ^ _p ^ 


( 11 ) 


with T'^(—p) := p),'0(p)) and 


T(p) ; = 


V'(p) 


The “Tr” denotes sum over momenta and internal indices. 
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FIG. 1: The quantum corrections to the four-Fermi interactions. The diagram surrounded by the 
dashed line is the large-leading term. The arrow denote the direction of particle number flow. 


Let us briefly explain how to introduce the RG flow equation of the four-Fermi cou¬ 
pling constant. The modihed inverse propagator -|- i?A can be divided into the held- 
independent part and the held-dependent part, 

+Ra = Sa + Va[% (12) 


with 


and 



Va: = 


7 

5TT 


G 


2N 
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of C ^ gT\ 

^ + S V j ’ 


where we employ the following notations 


=|l + z/i7o + ;irA(p) = {po + i^I,p{l + rAip)))^lf,, 
i>~ =^- ipio + i)rA{p) = (po - ip, P(1 + rA{p)))f, 7/., 
^ + ^7*y(^^75^)}, 

B = + (f75'0)i(^^7^)i}, 

c = 

V = {ipi'ipj + (275-0) i(i7®'0)j}. 


(13) 


(14) 


(15) 

(16) 

(17) 

(18) 
(19) 


In order to obtain the RG flow equation, we power expand the Wetterich equation with 
respect to Va[T] as follows, 

a^FA = -TT[dt ln(5A + Va[T])] 

= -Tr[4 ln(5A)] + Tr[4(5^'VA[^])] - ^Tr[4(5;^'VA[^])'] + • • • , (20) 


6 








FIG. 2; The threshold functions Iq (left) and Ii (right) on T — fl plane. 


where dt acts only on the A dependence of function R/^. We evaluate the third term as 
follow: 


-5Tr[4(5^'VA['I'l)=l = i 


/ d, 

r 1 1 1 

1 (27r)4 ^ 

[{p^y {p-y\ 


h! 

iV2 


(tr[7^Al7^Al]) + ■ • • , (21) 


where tr [7^717;,Al] = (tr[l spinor] ^ tr[lflavor])C)4Fermi = 4:N + (V>z75r/>)2}. This term 
yields factor N, thus it is the large-leading term. Other terms including the operators B, 
C and R do not yield it. For example, 


tr[7^T>7jC] = + (V'W5)KV’W^)m}] 

= - (^757m^)^}- (22) 

Using the Fierz transformation, we obtain the opertator {'ipipY + ('0*75'?/’)^. Although the 
nonleading terms yields not only the scalar-type operator but also other four-Fermi operators 
such as a vector-type operator ('07/A'0)^+(V’757/it/’)^, these terms are dropped in our truncated 
theory space. The possible ways of contraction between the operators are exhibited by the 
Feynman diagrams as shown in Fig. [1] In particular , the hrst diagram in the dashed box is 
the large-leading term. The large-A^ leading approximation neglects other diagrams than 
the hrst one. 

We obtain the RG how equation of the dimensionless four-Fermi coupling constant, 

dt9 = -2^ + ^{4/o(T, fi) - ^h{T, /i)}/, (23) 
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where the rescaled temperature T = T/A and rescaled chemical potential fi = p/A, obeys 
the following equations respectively, 


dtf = T, 

(24) 

dtp. = p. 

(25) 


The “form factor” Iq and Ji are the threshold functions of T and jl, which are given by 
the shell mode momentum integration of the large-leading term and the third diagram in 
Fig. dl respectively, and they are defined in Appendix The second and fourth diagrams in 
Fig.[T]do not contribute to the RG flow equation of g in the truncated effective action ([9]). We 


adopted the 3d optimized cut-off function 


22 


23l | as the regulator function R\ which makes 


it easy to integrate in momentum space with the Matsubara summation (see Appendix. ED- 
The large-leading approximated equation can be obtained by taking the limit A^ —)■ oo, 
which makes R term vanish. 

Let us now discuss the RG flows at hnite temperature and density in a qualitative manner. 
First, we discuss the hxed point structure of the RG equation. The non-trivial hxed point 
is given by 


9* =-^4-■ (26) 

hence, the hxed point moves with the RG evolutions of dimensionless temperature and 
density. When the threshold functions at hnite temperature and density are smaller than 
these with the vanishing temperature and density, we obtain > ( 74 = 0 ,/i=o- Thus, the 

chiral symmetry tends to be restored by thermal and hnite density ehects. 

Second, we briehy discuss the ehects of the large-A^ nonleading term at vanishing tem¬ 
perature and density. By taking the limit T —)■ 0 and p —)■ 0, we obtain 

9,g = -2(, + i (^4 - (27) 

where the factor 4 in the leading term is caused by the trace of the spinor space. Note 
that the diherence of coefficients, ^g"^ here compared with 2g'^ in Eq. (|5]), comes out of the 
diherence of cut-oh schemes. The nonleading term has negative sign, and it suppresses 
DySB. 

Finally, we discuss hnite temperature and density ehects. We show the threshold functions 
on T — fl plane in Fig. |2J The solutions of the RG equations of T and ft become simply the 








g g 


FIG. 3: The RG flows in case of the large-A^ leading at p, = 0.4 (left) and the large-A^ nonleading 
with N = 1 at p = 0.22 (right). 


exponentially-growing solutions with increasing t, thus, the ratio p/T is a scale independent 
constant, /r/T. Therefore the RG flows of T and p move on the straight line with the 
slope /i/T on T — p plane. Note that the threshold functions Jq and R with the vanishing 
temperature have the singular point at the Fermi surface A = /i: 


/o(0, /i) 

li{0,p) 


- 1 ), 

2(1 IaP ^ (T^ (5 ■ ^ 


(28) 

(29) 


where 6{p — l) is the step function. The function Jq at hnite dimensionless chemical potential 
and the vanishing temperature is constant and positive value. In particular, this function 
vanishes in the region A < /i. On the other hand, the function Ji depends on p in arbitrary 
scale. Furthermore, this function becomes quite large at the Fermi surface, therefore, the 
large-nonleading term is stronger than the leading term. This means that the RG flow of 
the four Fermi coupling constant with the large-nonleading term tends to go to positive 
inhnity at low temperature and high density. In the next section, we numerically evaluate 
the RG flow equations. 
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B. Numerical analysis 


We solve the RG flow equations of the inverse four-Fermi coupling constant g. The 
numerical results of the RG flows on g — fi plane are shown in Fig. |3l The left-hand side 
panel and the right-hand side panel are the large-A^ leading approximation case and the 
nonleading extended case, respectively. Here, we used the initial value: g = 0.3 where 
the RG flow at the vanishing temperature and density goes to the negative region in the IR 
limit. First, let us discuss the large-leading case. The RG flows of g at higher temperature 
cases (T = 0.2Ao and T = O.ISAq) go to positive infinity. The flow at T = O.IAq goes to 
negative region crossing the origin ^ = 0, and comes back to positive region. Eventually, 
this flow goes to positive infinity at the IR scale. The RG flows at lower temperature cases 
(T = O.OIAq and T = O.OOIAq) go to negative region at the IR limit. 

These flows can be evaluated only by the RG equation for the inverse four-Fermi coupling 
constant. If we solve the RG equation for the four-Fermi coupling constant Eq. fl25]) . these 
flows stop in the middle of the RG scale. Solving the RG equations of the inverse four- 
Fermi coupling constant is useful to discuss the flows after DySB, although it is not strictly 
authorized for large-nonleading case. It should be noted that to discuss RG flows towards 
positive infinity at IR scale do not always imply the restoration of the chiral symmetry. 
The RG flow of the inverse four-Fermi coupling constant corresponds to the curvature at 
origin of the effective potential with respect to the expectation value a ~ (V’V’)- case of 
the first order phase transition, the curvature at origin may remain positive, even after the 
DySB. We cannot distinguish between the restoration of chiral symmetry and the first order 
DySB. In other words, we see the effective potential only at the neighborhood of origin. The 
chiral phase transition including the first order phase transition should be investigated by 
searching the global minima of the effective potential by using the weak solution method or 
the re-bosonization method. Although our method can be applied only to the second order 
phase transition, our results may be used as references for advanced analysis of the chiral 
phase transition. 

Second, let us discuss the nonleading extended case. The RG flows at low tempera¬ 
ture (T = O.OIAq and T = O.OOIAq) go to positive infinity at the IR limit. These behaviors 
are quite different from the large-A^ leading case. The large-A^ nonleading term becomes 
larger than the leading term because the threshold function R has the singularity at /x = A. 
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FIG. 4: The chiral phase structure given by the large-A^ leading case (blue and green line) and 
nonleading case (red and pink line). We used go = 0.3 as the initial value and = 1 for the 
nonleading case. The phase boundaries assume the second order phase transition. 

The chiral phase diagram on /r/Ao — T/Aq plane is shown in Fig. HJ We show the phase 
boundaries evaluated by the blowup solution of Eq. fl2^ and the RG equation of the inverse 
four-Fermi coupling constant. Here we call the solution of the RG equation for the inverse 
four-Fermi coupling constant simply the inverse solution. 

First, we discuss the difference of phase boundaries between the blowup solution case 
and the inverse solution case. Since the blowup solution stops the RG flow in the midst 
of scale, we cannot evaluate it after the chiral symmetry breaking. Then we notice that 
the phase boundary with the inverse solution is correct. In high temeperature and low 
density region, both boundaries become same. By contrast, in low temperature and high 
density region, the phase boundaris evaluated by the blowup solutions become larger than 
the inverse solutions. That is, it is especially important to evaluate the RG equation of the 
inverse four-Fermi coupling constant in this region. 

Next, we discuss the difference of phase boundaries with the inverse solution between 
the large-leading case and nonleading case. The Dy^SB region becomes smaller in case of 
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the large-iV nonleading. We find drastic difference of phase boundaries at low temperature 
and high density region since the large-iV nonleading effect strongly suppresses the DySB 
as we previously mentioned. The chiral susceptibility, the curvature at the origin of mesonic 
potential ~ is sensitive to the nonleading effects at low temperature and high 

density. 

Note that the phase boundaries shown in Fig. 0] denote second order phase transition. 
As mentioned below Eq. ([H]), the inverse four-Fermi coupling constant corresponds to the 
curvature of the effective potential at the origin and its RG equation describes the change 
of the curvature. In the present work, the phase boundary is determined by the sign of the 
inverse four-Fermi coupling constant in IR limit. The broken phase potential with double 
well form has the negative curvature of the effective potential at the origin. The symmetric 
phase has a single well form and the positive curvature of the effective potential at the origin. 
Thus, the change from the negative curvature to the positive curvature mean the second 
order phase transition at least near the origin. One of important statements in our study 
is that these pictures can be described by the RG flow of the inverse four-Fermi coupling 
constant. 

On the other hand, when the first order phase transition occurs, there might be no global 
minimum at the origin of the effective potential even in case that the curvature of the 
effective potential at the origin changes to be positive. Therefore, it does not always mean 
we evaluate the global minimum of the effective potential in the present work. That is, there 
is possibility of occurring the first order phase transition at larger density region than the 
second order boundary evaluated in present paper. 

IV. SUMMARY AND DISCUSSION 

In this work, we have studied the RG flows of the four-Fermi coupling constant in NJL 
model at finite temperature and density by using the FRG. We solved the RG equations of 
the inverse four-Fermi coupling constant. This manner allows to integrate the equation up 
to the IR scale without the difficulties of the divergence. We have investigated the impact 
of the large-A nonleading term on the RG flow equation. The large-iV nonleading term has 
the singularity at the Fermi surface scale A = /i with the vanishing temperature, i.e. the 
quantum fluctuation due to the large-iV nonleading term becomes quite large, which makes 
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FIG. 5: The quantum corrections to the four-Fermi interactions generated by the Yukawa inter¬ 
action. The crossed ladder diagram (right) at the zero temperature and finite density has the 
singularity at the Fermi surface. 

the system more symmetric, as seen in Fig. H] 

In this analysis, the second order phase transition is investigated. In order to investigate 
the hrst order phase transition with the critical end point in the NJL model, the global 
minima of the effective potential with respect to the chiral order parameter must be studied 
by using more technical methods such as the re-bosonization method or the weak solution 
method. Currently, since the weak solution method can be applied to only the Fermionic 
large-approximated case, we need to employ the re-bosonization method. For example, a 
truncated effective action of the NJL model with bosonization is given by 

r 

Fa0 ] = y dr j +/i7o)'0 + ^a'0(o-+ • vrys)'!/’ 

+ + ^a(0^)| , (30) 

where 0* = (cr, tt). The four-Fermi interaction is generated by the Yukawa interaction 
in this action through the diagrams shown in Fig. O The contribution from the crossed 
ladder diagram at hnite density has the singularity at the Fermi surface. The impact of this 


singularity to the chiral phase diagram on T — /i plane should be investigated. The paper 


M 


25[ | which study QCD at vanishing temperature and density by using re-bosonization method 
indicates that the Yukawa coupling constant always gets the same value in the IR scale, 
i.e. the Yukawa coupling constant at the IR scale does not depend on the initial values. 
However, the impact of the singularity to the phase boundary is non-trivial. The analysis 
of the effective action flHUD will be presented elsewhere. 

On the other hand, the weak solution method is also powerful to analyze the dynamical 
chiral symmetry breaking 20|. This method is just developing to apply to many systems. 
Therefore, the present work becomes one of valuable benchmarks. We are going to clearly 
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discuss the relation between our results obtained in this paper and the phase diagram given 
by the weak solution method in future works. 
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Appendix A: Threshold functions 


In this appendix, we show the threshold functions which are obtained by the shell mode 
integration. We hrst introduce the cut-off prohle function Ra (p). In case of hnite tempera¬ 
ture and density, we employ the 3d optimized cut-off function j22|, l23| . 


-Ra(p) = *1^ 


A-d 


P\ ) 

=: ii) rA{p), 


9{1 - \p\/A) 


(Al) 


where 6{1 — \p\/A) is the Heaviside step function. 


e{i - \p\/A) = 


1 (|p|<A), 
0 (|p|>A). 


(A2) 


We dehne the following threshold functions. 


OO « 

k(T,,i-,A) = T J 


n=—oo 

OO 


d^p 

(27r)3"' [(po -K ip)2 + p‘^(l + rA(p))2] ’ 




(A3) 


n 

h(T,,i-,A) = T J 


d^p 

(27r)3"' [(po -K ip)2 + p2('x + rA(p))2][(po - ipY + p‘^{l + rA(p))2] ’ 




(Po + *h)(Po - ilA + + rMf 


(A4) 
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where dt = dtr^-dr andpo = (2n + l)7rT. We can analytically calculate the above integration 
and summation and we obtain 


Jo(T,/i;A) = 


A2 

T 


i-nA + Q-n-)+|;K + n-l 


o ;—>-1 


/i(T,/i;A) = — 


A2 ~ 

A^ r 1 /I \ 1/1 ^ 1 

F [(IT^ 


(AS) 


d 1 

n+ + 


9 


1 — jldu 


n_ 


A2 ~ 


(A6) 


Here we have used the formula, 


E 


1 

(2n + l) 27 r 2 + f] 2 El 


1 



(A7) 


where n± is the Fermi-Dirac distribution function: 


with E± = A± fi and (3 = 1/T. 


n± 


1 

exp(/3E±) + 1 


(A8) 
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